ADA  0  86825 


Technical  Report  535 


THE  DETECTABILITY  OF  A  MINE 
USING  REPLICA  CORRELATION  (RC) 
AND  ECHO-ECHO  CORRELATION  (EEC) 


LU 


NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CALIFORNIA  92152 


NOSC  TR  535 


NAVAL  OCEAN  SYSTEMS  CENTER,  SAN  DIEGO.  CA  92152 


A  N 


ACTIVITY  OF  THE  NAVAL  MATERIAL  COMMAND 
SL  GUILLE,  CAPT,  USN  HL  BLOOD 


Commander 


Technical  Director 


ADMINISTRATIVE  INFORMATION 


The  work  reported  herein  was  conducted  oyer  the  period  December  1979  to 
March  1980  in  support  of  CNO  Project  S0239-AS,  directed  by  NAVSEA  Code  63X1. 

The  author  wishes  to  express  his  appreciation  to  Mr.  S.  P.  Pitt  and  Mr. 

J.  Brady  of  ARL/UT  for  their  help  and  suggestions. 


Released  by 
RA  McLennan,  Head 
Sonar  Systems  Division 


Under  authority  of 
RD  Thuleen,  Head 

Weapon  Control  and  Sonar  Department 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PACE  (Mian  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 


NOSC  Technical  Report  535  (TR  535) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT’S  CATALOG  NUMBER 


4-r^TtrbE-(ond.Subtllto) - - - - - — . . .  .  . . 

t  Q>  J  The  Detectability  of  a  Mine  Using  Replica  Correlation  (RC)  and  I 

I  Echo-Echo  Correlation  (EEC)  |  * —  ""  ' 


».  CONTRACT  OR  GRANT  NUMBERfa; 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 


■  I.  CONTROLLING  OFFICE  NAME  AND  AODRESS 


Chief  of  Naval  Operations 
Washington,  DC  20350 


JJU  report  date 
//  Aprt»8gJ 


Washington,  DC  20350  ,3,  number  of  pages 

42 

<4.  MONITORING  AGENCY  NAME  ft  AODRESSf It  dlltarant  trom  Conttolllnt  Ottlce)  15.  SECURITY  CLASS,  (oi  thla  report) 


Unclassified 

15a.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (ol  thla  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abatract  entered  In  Block  30,  It  dlllarant  from  Rapotl) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Conllnua  on  ravaraa  alda  It  necaaeary  and  Identity  by  block  number) 


20.\^ABST R ACT  (Continue  on  ravataa  alda  It  nacaaaary  and  Identify  by  block  number) 

The  detectability  of  a  small  target  by  means  of  active  sonar  must  be  discussed  in  terms  of  waveform  characteristics 
and  the  signal  processing  used  to  extract  the  echo  from  interference.  This  report  presents  the  target  signature  of 
a  mine  for  large  bandwidth-time  product  (BT)  waveforms,  and  discusses  the  losses  suffered  when  the  echo  is 
processed  by  Replica  Correlation  (RC)  and  Echo-Echo  Correlation  (EEC).  It  is  shown  that  RC  “splitting” 
losses  increase  with  bandwidth,  while  EEC  losses  are  relatively  independent  of  bandwidth  and  are  much  lower 
than  RC  losses.  The  detectability  of  this  mine  for  a  given  scenario  of  false  alarm  rate,  number  of  beams,  sonar 


platform  motion,  and  processing  type  is  also  included.  , 

fS 


DD  |  JAN  ^73  1473  EDITION  OF  I  NOV  65  IS  OBSOLETE 


S/N  01 02-LF-01 4-6601 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (IWian  Dal*  Motored) 


jsoy 


UNCLASSIFIED 


SECURITY  CLAIMFICATION  OF  THIS  PAGE  <1Wi»n  Dttm  Bntirrd) 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGEOWnn  D*t»  Bnfr*d) 


CONTENTS 


I.  INTRODUCTION  . . .  page  3 

II.  TEST  PROCEDURE  ...  4 

III.  SIGNAL  PROCESSING  AND  DATA  REDUCTION  ...  6 

IV.  RESULTS  ...  28 

V.  DETECTION  OF  MINES  ...  33 

VI.  DISCUSSION ...  38 

Correlation  ...  38 
Target  strength ...  38 

VII.  CONCLUSIONS  ...  39 

VIII.  RECOMMENDATIONS  ...  39 

IX.  REFERENCES  ...  40 

X.  APPENDIX:  ECHO-ECHO  CORRELATION  ...  41 


1 


I.  INTRODUCTION 


This  paper  is  concerned  with  the  acoustic  detection  of  underwater  mines  from  the 
standpoint  of  signal  characteristics  and  signal  processing.  The  work  that  has  been  done  in 
the  field  has  relied  primarily  on  target  response  to  CW  pulses,  with  the  echo  processed  by  a 
simple  rectifier  and  averager  (Refs.  1,  2).  With  modem  microprocessors,  it  is  now  feasible  to 
employ  more  complex  signals  such  as  linear  frequency  modulated  (LFM)  pulses  and 
pseudo-random  noise  (PRN)  pulses  that  offer  greater  Doppler  tolerance,  better  noise  rejec¬ 
tion,  and  much  greater  signal  energy  than  is  possible  with  short  CW.  The  use  of  such 
signals,  however,  presents  its  own  problems,  particularly  as  regards  ocean  multipath  and 
target  structural  response  to  the  incident  acoustic  energy.  Since  adequate  data  were  not 
available  concerning  the  response  of  actual  mines  to  these  types  of  pulses,  it  was  decided 
to  conduct  an  experiment  to  determine  the  response  of  a  particular  mine  (Ref.  3)  to  several 
representative  waveforms.  The  test  was  devised  and  conducted  as  a  joint  project  between 
the  Applied  Research  Laboratory  (ARL)  at  the  University  of  Texas,  and  the  Naval  Ocean 
Systems  Center  (NOSC),  San  Diego,  California. 

One  of  the  measures  customarily  required  by  the  sonar  engineer  is  the  reflectivity 
index  of  a  targe ‘  generally  referred  to  as  the  “target  strength.”  Practically,  this  would  be  a 
ratio  of  the  received  echo  power  to  the  transmitted  power,  once  interfering  effects  such  as 
propagation  loss  and  absorption  are  accounted  for.  One  problem  that  arises,  the  one  perti¬ 
nent  to  this  discussion,  is  the  difficulty  in  making  quantitative  measurements  of  the  received 
echo.  As  shown  later,  the  echo  structure  is  highly  dependent  on  the  aspect  angle  of  the  tar¬ 
get,  and  as  such,  would  be  very  difficult  to  measure  in  an  operational  environment. 

The  detectability  of  an  acoustic  target  when  usirg  relatively  complex  waveforms 
is  enhanced  by  some  form  of  signal  processing.  Two  types  of  coherent  processing  are 
discussed  in  this  paper,  namely  Replica  Correlation  (RC)  (Refs.  4,  5)  and  Echo-Echo 
Correlation  (EEC)  (Ref.  5),  in  terms  of  their  ability  to  increase  the  detectability  of  this 


mine. 


II.  TEST  PROCEDURE 


The  Target  Strength  Measurement  Facility  at  ARL  (Ref.  6)  was  used  to  conduct  this 
experiment.  The  target  was  held  as  shown  in  figure  1 ,  with  the  major  axis  in  the  horizontal 
plane  of  the  hydrophones.  The  cuntrol  mechanism  allows  either  continuous  rotation  of  the 
body,  or  rotation  to  fixed  aspect  angles,  as  desired.  For  this  experiment  the  target  was 
rotated  to  twenty  different,  discrete  aspect  angles.  The  target  geometry  is  shown  in  figure  2, 
with  the  aspect  angles  investigated  listed  on  that  figure.  Note  that  the  “zero”  aspect  angle 
for  this  target  occurs  at  broadside,  a  peculiarity  of  the  AllL  rotator. 

The  transmit  pulses  were  generated  by  a  programmatic  device  designed  at  ARL.  The 
eight  pulse  codes  (PC’s)  used  are  listed  in  table  1.  At  the  conclusion  of  the  eightieth  trans¬ 
mission,  the  target  was  moved  to  the  next  aspect  angle,  and  the  sequence  repeated. 


PC 

B(kHz) 

T(ms) 

Type 

% 

l 

1 

100 

LFM  single 

2 

2 

100 

LFM  single 

3 

4 

100 

LFM  single 

4 

1 

3X24.6 

LFM  triplet 

5 

2 

3X24.6 

LFM  triplet 

6 

4 

3X24.6 

LFM  triplet 

7 

4 

!00 

PRN  single 

8 

4 

3X24.6 

PRN  triplet 

B  =  bandwidth,  T  =  pulse  length;  Center  Frequency  =  27 .6kHz 
Table  1.  Experimental  waveforms. 

In  addition  to  the  target,  a  calibrated  probe  was  placed  at  a  known  distance  from  the 
hydrophone  for  one  series  of  transmissions.  The  data  obtained  were  intended  for  measuring 
path  loss  and  target  strength  and  system  distortion. 

The  acoustic  returns  were  electronically  gated  so  that  a  small  portion  of  pre-echo 
noise,  the  echo,  and  some  trailing  noise  were  recorded.  The  gated  returns  were  sampled 
with  a  12-bit  digitizer  and  the  samples  stored  on  magnetic  tape. 

Finally,  a  series  of  pulses  was  sampled  directly  at  the  output  of  the  signal  generator, 
and  the  samples  were  stored  as  before  for  use  as  replicas  in  later  signal  processing. 
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Figure  2.  Target  geometry. 


III.  SIGNAL  PROCESSING  AND  DATA  REDUCTION 

The  pulses  discussed  in  this  paper  all  have  bandwidth  -  time  products  (BT)  much 
greater  than  a  CW  pulse.  The  use  of  these  pulses  in  acoustic  detection  is  enhanced  by  some 
form  of  processing  which  reduces  the  BT  product  to  a  lower  value.  This  may  be  accomplished 
by  incoherent  processing,  such  as  rectification  followed  by  averaging  or  by  coherent  process¬ 
ing,  such  as  correlation.  The  latter  approach  is  considered  here,  and  is  limited  to  Replica 
Correlation  and  Echo-Echo  Correlation.  RC  is  the  “optimum”  filter  for  an  echo  imbedded  in 
Gaussian,  white  noise,  when  that  echo  exactly  replicates  the  transmitted  pulse,  while  EEC  is 
a  relatively  new  process  which  tends  to  account  for  environmental  and  target  induced  distor¬ 
tions  of  the  echo.  A  summary  of  the  EEC  process  is  given  in  the  appendix.  In  the  RC  case, 
for  an  ideal  situation  of  echo  perfectly  matched  to  the  transmitted  pulse,  the  peak  of  the 
correlation  of  the  echo  should  represent  the  acoustic  power  contained  in  the  echo.  For  EEC, 
an  equivalent  result  should  hold,  regardless  of  the  mismatch  between  transmit  pulse  and  echo. 

The  data  reduction  falls  into  three  categories.  First,  ten  echoes  for  each  pulse  code 
at  each  aspect  angle  were  envelope  detected,  squared,  and  averaged  together.  The  results 
are  shown  in  figures  3-10.  The  average  power  at  each  aspect  angle  was  determined  as  the 
average  mean  square  amplitude  of  the  echo  set  prior  to  detection.  Second,  each  echo  was 
correlated  with  its  appropriate  replica.  The  resulting  curves  were  envelope  detected,  squared, 
and  averaged  together  as  was  done  for  the  echoes.  The  resulting  curves  are  shown  in  figures 
11-18.  Note  that  the  time  axis  has  been  expanded  to  show  30ms  about  the  main  peak.  Third, 
the  EEC  process  was  performed  for  the  appropriate  echoes  (PC  4,  5,  6,  8),  the  results  envelope 
detected,  squared,  and  averaged,  as  was  done  previously.  The  results  are  shown  in  figures  19-22. 
The  plots  present  the  same  expansion  of  the  peaks  as  do  the  RC  plots,  but  the  abscissa  must  be 
interpreted  as  discussed  in  the  appendix.  The  vertical  scale  for  the  three  sets  of  plots  is 
indicated  on  each  plot. 

Finally,  the  processing  gain  available  for  the  two  types  of  correlation  was  explored. 

To  that  end,  ten  independent,  Gaussian,  bandlimited  white  noise  populations,  each 
normalized  to  an  integrated  power  spe  tral  density  of  unity,  were  added  in  succession  to  the 
successive  echoes  for  the  desired  pulse  codes,  for  selected  aspect  angles.  The  input  signal  to 
noise  ratio  was  computed  as 

SNi  =  Power /o?  (1) 

where:  “Power”  is  taken  from  the  echo  plots,  figs  3-10  as  appropriate;  of  is  the  average 

variance  of  the  input  noise  populations  (a  multiplicative  scale  factor). 

The  output  signal  to  noise  ratios  were  computed  as 

SNo  =  Poak/oQ  (2) 

where:  “Peak”  is  the  highest  (average)  output  of  signal  plus  noise;  a~  is  the  average 

variance  of  the  noise  only  correlation  output,  prior  to  envelope  detection. 
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The  theoretical  SNo’s  (Ref.  5)  are,  for  RC 

SNo  =  2BT(SNi),  (3) 

and  for  EEC 

SNo  =  2BT(SNi)2.  (4) 


Echoes 


TIME 


Figure  9.  Pulse  Code  7. 
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Figure  10.  Pulse  Code  8. 


IV.  RESULTS 


The  most  striking  feature  of  the  echo  plots  (Figs.  3-10)  is  the  extreme  variability, 
both  in  echo  shape  and  echo  power.  It  should  be  recalled  that  the  power  in  an  echo  was 
calculated  as  the  mean  square  amplitude  (prior  to  envelope  detection).  For  the  shapes 
shown  here,  it  is  likely  that  such  a  measure  is  not  optimum,  yet  it  remains  a  reasonable 
measure  against  which  to  compare  later  results.  In  an  operational  setting,  however, 
computing  the  power  in  this  way  first  requires  that  the  echo  be  identified,  and  the  echoes 
as  shown  might  prove  difficult  to  recognize  in  background  interference.  In  such  a  case,  some 
form  of  signal  processing  is  needed  to  aid  in  echo  identification. 

The  RC  plots  are  shown  in  Figs.  1 1-19.  Again  one  of  the  main  features  is  the 
extreme  variability  in  both  shape  and  peak  level.  For  an  ideal  echo,  that  is,  one  which 
matches  the  transmitted  pulse,  there  would  only  be  a  single  “spike”  in  the  RC  output  when 
the  echo  had  passed  through  the  processor.  A  good  example  is  shown  in  all  RC  plots  at 
an  aspect  angle  of  200  degrees  (compare  with  the  echo  plots  for  the  same  angle).  It  is  seen 
here  that  the  transmitted  pulse  has  essentially  generated  an  echo  from  a  single  reflector  on 
the  target.  In  the  case  of  PC4  (Figs.  6  and  14),  the  echo  power  is  -27  dB  and  the  RC  peak  is 
-28  dB,  a  very  close  agreement.  At  the  other  extreme,  however,  consider  the  results  at  an 
aspect  angle  of  1 10°.  Again  for  PC4,  the  echo  power  is  -38  dB,  while  the  RC  peak  is 
-43  dB,  a  difference  of  5  dB.  For  the  higher  resolution  pulses,  the  differences  tend  to  be 
somewhat  greater. 

It  was  anticipated  that  the  higher  resolution  pulses  would  also  tend  to  reveal  more  of 
the  target  structure  than  the  smaller  bandwidth  pulses.  In  fact,  this  generally  was  not  true. 
PCI  reveals  very  nearly  as  much  structural  detail  as  does  PC3,  except  at  aspect  angles  near 
100°.  It  is  thus  concluded  that  for  this  target,  spatial  resolution  less  than  1  yard  (B~l  1Hz) 
is  not  necessary  when  using  RC  processing.  However,  there  may  be  a  requirement  for 
larger  B  to  enhance  the  processing  gain  against  reverberation. 

The  PRN  echoes  and  RC  curves  (PC  7,  8)  are  quite  disappointing  in  two  respects. 
First,  these  echoes  had  a  much  poorer  signal-to-noise  ratio  (not  shown  in  the  figures)  than 
the  other  PC’s,  caused  by  the  necessity  of  transmitting  a  reduced  signal  power  in  order  to 
accommodate  large  signal  amplitudes  without  distortion.  Second,  and  perhaps  because  of 
the  first  effect,  the  RC  plots  are  very  ragged  and  lack  the  abrupt  dropoff  on  either  side  of 
the  main  response  lobe,  as  is  demonstrated  by  the  LFM  RC  curves.  It  would  thus  appear 
that  this  type  of  PRN  pulse  followed  by  RC  would  not  be  a  suitable  waveform  in  a  noise- 
limited  environment. 

The  EEC  plots  are  shown  in  Figs.  19-22.  A  major  difference  between  these  and 
the  RC  curves  is  the  single  peak  located  precisely  at  the  center  of  the  EEC  segment.  Note 
also  that  the  curves  are  symmetrical  within  a  segment.  Unlike  the  RC  plots,  which  showed 
little  increase  in  structural  detail  with  increasing  B,  the  EEC  plots  show  definite 
improvement  in  this  regard. 

Of  particular  interest  is  the  EEC  result  for  the  PRN  pulses,  PC  8  (Fig.  22).  Compare 
this  plot  with  the  EEC  result  for  PC  6  (Fig.  21)  and  with  the  RC  result  for  PC  7  (Fig.  17). 

It  is  seen  that  the  EEC  process  produces  results  which  are  very  nearly  identical  for  both 
LFM  and  PRN,  and  which  are  considerably  more  usable  than  RC  for  PRN  pulses. 
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An  analysis  of  Figs.  23-27  reveals  that  the  loss  between  power  and  peak  for  RC 
increases  with  bandwidth,  while  the  EEC  values  are  relatively  unaffected.  This  result  is 
emphasized  in  Figs.  28  and  29,  which  show  a  comparison  for  RC  and  EEC,  respectively, 
between  an  average  (over  PC)  power  level  and  the  correlated  peaks.  On  these  figures,  the 
loss  shown  at  each  aspect  angle  reflects  that  loss  which  occurred  between  power  and  peak  for 
the  same  PC,  as  shown  in  the  previous  five  figures.  These  results  indicate  that  EEC  is  able  to 
recombine  the  energy  distributed  throughout  the  echo,  wliile  RC  is  less  able  to  do  so. 

The  probe  “echoes”  were  correlated  with  the  stored  replicas,  and  the  results 
compared  with  the  replica  autocorrelations.  No  distortion  or  losses  of  any  kind  were 
discernible,  which  indicates  that  the  transmit  system  did  not  distort  the  waveforms. 

It  was  hypothesized  that,  in  an  ideal  (Gaussian,  white)  noise  background,  any 
loss  in  the  (average)  correlated  output  would  be  due  solely  to  target  structural  influences, 
even  with  a  spatially  stationary  target.  This  was  tested  in  the  previously  mentioned  signal 
to  noise  analysis  by  using  as  the  input  power  value  the  value  of  the  respective  peak  (RC  or 
EEC)  shown  in  Figs.  23-27,  rather  than  the  power  value  shown  there.  It  was  found  that 
the  measured  results  then  were  in  close  agreement  with  the  predicted  values  from  Eq. 

(3)  or  (4),  as  appropriate.  It  is  thus  concluded  that  the  hypothesis  is  correct,  and  that  the 
theoretical  processing  capability  of  either  RC  or  EEC  is  valid  provided  the  input  signal 
level  can  be  determined. 

Finally,  it  should  be  noted  that  Figs.  23-29  represent  the  target  strength  of  this  mine 
for  the  aspect  angles  indicated.  It  must  be  emphasized  that,  in  an  operational  environment, 
the  raw  echo  could  not  be  detected  except  at  very  high  echo  to  interference  levels,  so  that 
the  sonar  operator  would  have  to  rely  on  some  form  of  echo  processing  to  provide  better 
detectability. 


ASPECT  ANGLE,  deg. 


Figure  23.  Comparison  of  Echo  Power,  RC  Peak  and  EEC  Peak 
for  Pulse  Code  4. 


*based  on  calibrations  supplied  by  ARL/UT. 
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Figure  26.  Comparison  of  Echo  Power,  RC  Peak  and  EEC  Peak 
for  Pulse  Code  8. 
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Figure  27.  Comparison  of  Echo  Power  for  PC’s  4, 5, 6, 8. 
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DIFFERENCES  IN  ECHO  POWER  AND  RC  PEAK  AS  A 
FUNCTION  OF  ASPECT  ANGLE  AND  BANDWIDTH 
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Figure  28.  Differences  in  Echo  Power  and  RC  Peak  as  a  function  of 
aspect  angle  and  bandwidth. 
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Figure  29.  Differences  in  Echo  Power  and  EEC  Peak  as  a  function  of 
aspect  angle  and  bandwidth. 


V.  DETECTION  OF  MINES 


The  purpose  of  this  section  is  to  discuss  the  detection  of  mines  under  an  assumed 
scenario,  and  then  relate  the  discussion  to  the  results  of  the  previous  section. 

Assume  the  following  sonar  requirements:  resolution  of  1  yard;  pulse  length 
t  —  75  ms  for  RC  processing,  r  =  25  ms  for  EEC  processing;  N  =  20  beams  in  azimuth;  own 
velocity  uncertainty  is  ±1  kt;  bearing  uncertainty  is  ±1  beam.  The  false  alarm  problem  is 
somewhat  unusual  and  is  as  follows:  there  is  a  50%  probability  of  no  false  alarms  in  1 
hour.  The  observation  time,  T,  is  actually  modified  by  limiting  the  attention  of  the  observer 
to  the  last  2/3  of  each  ping  history.  Further,  the  decision  rule  will  be  based  on  two  out  of  three 
occurrences,  but  this  will  be  modified  as  indicated  later. 

The  resolution  capability  requires  a  bandwidth  of  about  B  =  1  kHz.  Thus  2Bt  =  50 
for  EEC,  and  2Br  =  150  for  RC. 

Figure  30  is  an  idealized  plot  of  the  probability  density  functions  (pdfs)  of  the 
noise-only  and  signal-plus-noise  distributions.  The  probability  of  false  alarms  is  p',  so  that 
p  =  1  -  p'.  Now 

prob  [  1  or  more  out  of  k  events]  = 

1  -  prob  [zero  out  of  k]  = 

l--[(l-p')kL  (5) 

which  is  the  binomial  distribution  for  zero  occurrences  out  of  k  trials.  Define  k  to  be 


k  =  N  •  B  •  T  •  2/3  (6) 


The  problem  may  be  stated  as  follows:  If  there  is  a  50-percent  probability  of  one  or 
more  false  alarms  in  k  trials,  what  is  the  required  p'  of  a  single  trial? 

Thus 

0.5  =  1  -  [(l-p')k]  (7) 

or 

p'  =  1  -  (0.5)  ^  (8) 
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It  is  probable  that  the  decision  will  not  be  made  on  a  single  threshold  crossing,  but 
rather  on  m  crossings  out  of  k  possibilities.  This  again  is  a  binomial  distribution  problem,  as 
follows:  The  probability  of  m  occurrences  of  an  event  in  k  independent  trials  is  given  by 

Pw.m  "Ecfp" qk'm  (9) 


where 


q  =  1  -p 


rm  .. 
K 


k! 

ml  (k-m)! 


Assuming  two  out  of  three  trials  defines  the  decision  rule, 

Pk;m  =  2!  (3-2)1  ^  ^  ^  +  3!  (3-3)!  ^  ^ 

-  3(p')2 
=  p'  (Eq.  (8)). 


where  p'  is  the  required  probability  of  false  alarms  for  one  event  necessary  to  meet  the 
decision  rule. 

The  rate  of  exceeding  a  threshold  is  the  False  Alarm  Rate,  K.  When  B  is  large  an 
approximation  is  (Ref.  7) 

K'  =p'B.  (11 

Assuming  that  the  operator  is  observing  the  detected  envelope  of  the  correlation,  one  may 
use  the  following:  define  the  recognition  differential  to  be  (Ref.  5) 

d  =  201og[(D-aVff/2)/0.655or] 

D  =  [-2a2  2n(K'/B)]  X>2 
(B  in  Hz) 

(assume  a  =1).  ty 


Using  curves  from  Ref.  5,  one  obtains  from  d  the  50%  Signal  Differential  (SD)  or  signal-to- 
noise  ratio  (in  dB)  into  the  correlator  necessary  to  achieve  the  probability  of  false  alarms, 
p' ,  for  a  50%  probability  of  detection.  However,  for  EEC,  it  is  possible  to  reduce  K' 
simply  by  not  computing  the  correlations  except  about  the  zero  lag  position.  The  minimum 
K  for  EEC  is  about  K'/2B r.  When  the  target  or  own  ship  has  a  velocity  uncertainty,  that 
uncertainty  will  define  the  minimum  K  for  EEC.  For  the  stated  uncertainty  of  v  =  1  kt, 
one  has 
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AR  =vr  =  1.68  ft/sec  •  0.025  sec  =  4.2  X  10-2  ft 


(13) 


so  that 


5  =  AR/C  =  8.4  X  10"6  sec. 

EEC  compensates  for  this  time  uncertainty  merely  by  searching  the  correlation  output  for 
±5  seconds  about  the  zero  lag  position.  The  maximum  reduction  in  the  output  duration 
compared  with  the  input  duration  is 

a  =  28/r, 

however,  for  security,  let  us  search  56,  so  that 

a  =  1.68  X  10"3.  (14) 

Then  for  EEC 

k  =  K  'la.  (15) 

Using  the  above-mentioned  curves,  one  has  for  RC 


SD  =  -9.8  dB  (16) 

while  for  EEC 

SD  =  -5.3  dB  (17) 

At  this  point,  one  must  more  carefully  account  for  the  uncertainty  in  own  ship’s  velocity. 
The  assumption  of  ±1  knot  potential  error,  combined  with  an  illumination  of  a  given  sea 
volume  every  10  seconds,  causes  a  potential  range  uncertainty  of  ±6  yards.  Replica  correla¬ 
tion  provides  a  much  more  precise  estimate  of  target  range  than  docs  EEC,  and  is 
correspondingly  more  sensitive  to  uncertainties  in  velocity.  Since  ±6  yards  is  approximately 
equal  to  ±6  resolution  cells,  the  following  rule  is  assumed  for  RC:  Whenever  the 
from  any  one  resolution  cell  exceeds  a  threshold  in  one  transmission,  the  next  transmission 
will  require  the  observer  to  consider  6  resolution  cells  on  either  side  of  the  original  cell  as 
potentially  containing  the  echo  from  the  same  original  target.  On  the  third  transmission, 
the  observer  will  then  be  required  to  observe  a  total  of  25  cells.  The  situation  is  further 
compounded  by  the  beam  uncertainty.  Again  for  RC,  the  3-beam  rule  requires  that  39 
cells  be  searched  on  the  second  transmission,  while  125  cells  are  searched  on  the  third. 

See  figure  3 1 .  For  EEC,  on  the  other  hand,  the  inherent  range  uncertainty  is  cr  =  ±41  yards, 
so  that  the  velocity  uncertainty  is  transparent  to  an  observer  at  a  fixed  range.  However, 
two  transmissions  will  result  in  a  compounded  potential  error  of  ±12  yards,  which  is  29% 
of  the  inherent  uncertainty.  Thus  it  will  be  assumed  that,  for  EEC,  the  following  rule 
applies:  When  an  echo  is  identified  on  one  transmission,  the  second  transmission  will 
necessitate  the  observation  of  3  cells  (1  range  cell  by  3  beam  cells),  while  the  third 
transmission  will  require  the  observation  of  a  total  of  10  cells. 


Figure  31.  Multiping  range  and  bearing  uncertainty. 


The  probabilities  previously  calculated  did  not  include  the  velocity  uncertainties 
Thlw  restating  the  problem:  Given  a  probability  f)  of  a  single  occurrence  on  the  first  pn  g, 
what  is  the  probability  of  an  occurrence  in  the  second  ping  or  in  the  third  ping.  I  shou 
be  noted  that  an  occurrence  on  the  first  ping  presupposes  a  single  target  so  that  the 

mor“  ha„  one  cell  being  “occupied”  on  the  second  or  thnrd  pmg  .s  very 

small.  From  Eq.  (9),  for  RC, 


~3 9  p(l-38$) 

(For  the  second  ping) 


and 

p125;l=  125  *  pCl-p)124 
~  125  •  p( 1— 124p) 

(For  the  second  ping) 

So,  for  RC,  the  total  probability,  which  is  equal  to  p'  [Eq.  (8)],  is 
P'  =  P’P39;1  +  P*Pl25;0'Pl25;l 
39p2  +  125p2 
=  164p2 . 


(19) 


(20) 
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For  EEC,  we  have 

p3;1  =  3  •  p(l-p)2 
3p 

(For  the  second  ping)  (2  I) 

and 

P  10;1  “  10 'P 

(For  the  tliird  ping)  (22) 

Thus  the  total  probability  for  EEC  is 
p'  =  3p2  +  10p2 

=  13p2.  (23) 

The  results  may  finally  be  given  as  in  Table  2. 


_P_ 

(K '/«) 

d 

_SD_ 

RC 

9.37  X  50-6 

9.37  X  10~3 

14.7 

-8.7 

EEC 

3.33  X  10"5 

19.82 

7.5 

-4.5 

(Note:  a  = 

1.0  for  RC) 

Table  2.  Signal  differentials. 


These  results  indicate  that  EEC  requires  approximately  4  dB  greater  SNi  than  does  RC, 
for  a  perfect,  matched  echo.  However,  from  Fig.  23,  it  is  seen  that  RC,  on  the  average, 
suffers  3  ±  1.5  dB  “splitting”  loss  for  PC  4,  whereas  EEC  suffers  an  average  reduction  in 
power  recombination  of  1 .4  ±  0.5  dB.  However,  the  maximum  EEC  loss  for  this  pulse  code 
was  2  dB,  whereas  that  for  RC  was  5  dB.  In  addition,  the  difficulty  in  displaying  27,000 
output  samples  to  the  operator  will  probably  require  a  reduction  in  the  number  of  points 
through  “peak  picking”  or  “or-ing”  (Ref  8).  Such  a  reduction,  for  the  recognition 
differential  (d)  listed  in  Table  2,  will  entail  an  additional  loss  in  S.D.  for  RC  of  at  least 
several  dB.  This  would  place  the  SD  for  RC  about  equal  to  that  required  for  EEC. 


37 


VI.  DISCUSSION 


CORRELATION 

For  a  given  pulse  bandwidth,  EEC  requires  three  times  the  pulse  length  necessary  for 
RC  for  a  given  BT  product.  For  systems  which  require  short  transmit  times,  either  for 
electronic  or  security  reasons,  this  may  be  a  serious  difficulty.  EEC  does  not  provide  precise 
range  information  as  does  RC.  There  is  an  ambiguity  in  the  range  comparable  to  the  travel 
time  of  one  third  of  the  transmitted  triplet  pulse.  However,  the  resolution  capability  of 
EEC  is  the  same  as  for  RC.  EEC  has  the  following  advantages  over  RC: 

1 .  Simplicity  of  Processing.  Replica  correlation,  even  when  done  by  “fast  cor¬ 
relation,”  is  a  time  consuming,  and/or  hardware  consuming  operation.  The  array  pro¬ 
cessor  necessary  to  satisfy  system  requirements  is  likely  to  be  large  and  expensive.  EEC, 
on  the  other  hand,  may  be  accomplished  by  simple  time-domain  correlation  in  real-time, 
since  only  those  correlation  lags  necessary  to  account  for  unknown  Doppler  must  be 
computed. 

2.  Reduced  Display  Requirements.  It  is  estimated  that,  using  a  1-kHz  B  pulse  in 
an  operational  setting,  RC  processing  will  require  the  observation  of  27,000  spatial  reso¬ 
lution  cells  per  transmission  for  20  beams.  EEC  will  require  the  observation  of  only  900 
such  cells. 

3.  Energy  Recombination.  The  structural  properties  of  the  target  will  produce  a 

‘  split”  echo.  RC  is  not  capable  of  recombining  this  energy,  whereas  EEC  does  so  very  well. 
In  addition,  RC  shows  increased  decorrelation  with  increasing  bandwidth,  whereas  EEC  does 
not  seem  to  be  so  affected. 

4.  False  Alarm  Rate.  The  false  alarm  rate  for  EEC  may  be  greatly  reduced  simply 
by  not  computing  the  correlation  function  at  times  far  removed  from  the  zero  lag  (segment 
center).  This  has  the  effect  of  reducing  required  detection  thresholds,  thus  improving 
target  detectability. 

5.  Target  Motion.  For  the  waveforms  studied  here,  RC  processing  is  extremely 
sensitive  to  target  and/or  own-ship  velocity,  and  may  show  severe  distortion  in  the  presence 
of  uncorrected  accelerations.  A  more  appropriate  waveform  to  use  in  such  a  case  is  linear 
period  modulation  (LPM),  which,  for  RC,  will  provide  Doppler  tolerance.  However, 
accelerations  will  still  prove  highly  distorting.  For  EEC,  using  LPM  pulses,  the  results  will 
be  insensitive  to  Doppler,  and  highly  tolerant  of  constant  or  slowly  varying  accelerations. 

6.  Pseudorandom  Noise.  It  appears  that  PRN  pulses  of  the  type  used  here  will 
not  be  appropriate  for  RC  processing  in  a  noise  environment,  but  will,  aside  from  the 
reduced  transmit  power  mentioned  previously,  work  as  well  as  LFM  pulses  for  EEC 
processing. 

TARGET  STRENGTH 

With  high-resolution  pulses,  the  concept  of  target  strength  is  somewhat  ambiguous. 
One  must  discuss  this  measure  in  terms  of  the  signal  processing  to  be  used,  and  it  must  be 
clear  if  the  measure  is  an  average  over  some  spatial  extent  of  the  echo  or  is  a  peak  value. 
Assuming  that  coherent  echo  processing  is  operationally  required,  the  processor  must  decide 
whether  to  average  the  correlation  output  over  some  time  greater  than  the  basic  resolution 
of  the  pulse,  or  not.  It  is  the  author’s  contention  that  a  proper  averaging  time  is  totally 
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dependent  on  the  aspect  angle  of  the  target,  and  that,  without  such  information,  post¬ 
correlation  averaging  will  be  counterproductive.  However,  any  post  correlator/detector 
averaging  suitable  for  RC  would  be  identically  suitable  for  EEC,  so  the  comparison  of  the 
two  is  valid,  independent  of  this  additional  processing.  The  results  presented  here  reflect 
that  which  a  sonar  operator  would  encounter,  provided  the  background  interference  is  a 
noise  limiting  environment.  However,  for  large  bandwidth  signals,  it  is  felt  that  the  same 
results  will  apply  for  a  reverberation  interference. 

VII.  CONCLUSIONS 

An  investigation  of  tiie  acoustic  response  of  a  particular  mine  has  been  presented  in 
terms  of  aspect  angle  dependency  and  signal  processing.  It  has  been  demonstrated  that  EEC 
is  superior  to  RC  for  recombining  “split”  energy  from  this  mine  and  the  results  indicate  that 
EEC  would  be  operationally  superior  to  RC  in  mine  detection  in  a  noise-limited  environment, 

In  addition  to  investigating  this  one  target,  it  has  been  demonstrated  that  an  effec¬ 
tive,  rapid  methodology  is  available  for  obtaining  the  acoustic  signature  of  any  small  mine 
or  torpedo-like  object  for  any  pulse  type  conceivable.  Signal  processing  techniques  such  as 
simple  averaging,  RC,  and  EEC  are  available,  and  others  can  be  incorporated. 

VIII.  RECOMMENDATIONS 

A  similar  series  of  tests  should  be  conducted  in  reverberation  environments  to  deter¬ 
mine  the  relative  effect  on  RC  and  EEC. 
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X.  APPENDIX:  ECHO-ECHO  CORRELATION 

The  purpose  of  this  appendix  is  to  provide  a  brief  summary  of  the  EEC  technique. 
For  more  detailed  information,  the  reader  is  directed  to  Ref.  5. 

Consider  the  envelope  of  an  acoustic  transmission  such  as  is  shown  in  Fig.  Al. 

Note  that  the  individual  length,  r,  of  each  pulse  in  the  triplet  is  <  L.  Assume  that  the  range 
of  the  target  is  unknown.  Immediately  upon  completion  of  the  transmission  fat  time  3D. 
begin  partitioning  the  received  time  history  (RTH)  into  segments,  each  L  seconds  in  length. 
After  some  time  commensurate  with  the  range  of  the  target,  an  echo  set  such  as  shown  in 
figure  A2  will  be  received.  Note  the  arbitrary  alignment  of  the  echoes  vis-a-vis  the  segmen¬ 
tation.  At  any  time  prior  to  the  arrival  of  the  echo  set,  take  segment  X^-,  sample  it  and 
correlate  it  with  the  sampled  version  of  segment  Xj^  +  j,  as  indicated  in  Eq.  (Al). 

+L 

CK(X)  =  ^  X  XK <j+X)  XK+I<J>  •  (Al) 

j~L 

Next,  discard  segment  Xj^  and  correlate  Xj^+]  with  Xj£+2,  obtaining  Cjr+j(X).  Note 
that  if  segments  Xj^  and  Xj^+1  contain  only  noise,  then  the  expected  value  of  Cj^(X)  is 
zero.  The  same  is  true  if  X^+i  is  noise  only,  while  Xj^+2  is  signal  plus  noise.  However, 
if  Xj£+2  and  X^+3  both  contain  echo,  then  the  expected  value  of  Cj^fX)  is  the  auto¬ 
correlation  of  the  echo,  contaminated  with  noise. 

An  inspection  of  this  process  reveals  two  important  items: 

1 .  For  a  spatially  stationary  target,  the  peak  of  the  echo  correlation  falls  precisely 
at  Cj£+2  (0),  ie,  at  the  center  of  the  segment. 

2.  Based  on  item  1 ,  above,  it  is  only  necessary  to  compute  C£+2(X)  at  X  =  0, 
unless  the  target  is  moving.  In  this  case,  one  need  compute  Cj^fX)  only  for  values  of 
8  <  |  Xmax  |  =  Vr/c,  where  V  =  target  velocity  (relative),  C  =  sonic  speed.  Thus  the 
opportunity  for  a  false  alarm  is  reduced  in  relation  to  V. 

The  EEC  plots  shown  in  the  text  represent  I  Xmax  I  =  2.5  ms.  While  each  time  mark 
on  the  abscissa  represents  25  ms  of  elapsed  time,  the  resolution  between  time  marks  is  only 
5  ms.  Thus,  while  the  actual  elapsed  time  to  the  target  is  known  only  to  within  25  ms, 
the  basic  resolution  of  the  pulse  is  retained,  and  is  the  same  as  for  RC. 
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